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We present a detailed analysis of the ion hopping dynamics and the related nearby oxygen dynamics 
in a lithium meta silicate glass via molecular dynamics simulation. For this purpose we have 
developed numerical techniques to identify ion hops and to sample and average dynamic information 
■ of the particles involved. This leads to an instructive insight into the microscopic interplay of ions 

and network. It turns out that the cooperative dynamics of lithium and oxygen can be characterized 
^ ' as a sliding door mechanism. It is rationalized why the local network fluctuations are of utmost 
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importance for the lithium dynamics. 
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I. INTRODUCTION 



As a well known characteristic feature of ion conducting glasses the dynamics of ions 



and network atoms are strongly_ 
glass transition temperature 



decoupled at low temperatures and in particular below the 



. Although by definition for temperatures T < Tg the 
network dynamics is mainly limited to localized fluctuations, molecular dynamics simulations 
of alkali silicate systems with an immobilized network give clear evidence that the network 
d.„a.,. a ..nmc. „.pac. on .Ke .on d_ Q. A ...Ua. beKa™. 

has already been predicted by Anderson and Stuart who have assumed that the opening of 
diffusion pathways between adjacent ionic sites through the silicate matrix supplies a major 
contribution to the activation energy of the ion diffusion in alkali silicate glasses 0, 0]. In 
recent work it has been proposed that at least for mixed-alkali systems a strong coupling 
of alkali and network dynamics should be present to rationalize experimental findings jlol |. 
Furthermore for models of ion dynamics it is particularly important to know the nature of 
the network relaxation if an ion is leaving its site 

In [21 we have recently shown that for a lithium silicate system with an artificially im- 
mobilized silicate network the lithium dynamics changes in many aspects compared to a 
system with unchanged network mobility. A dramatic slowing down of the ion dynamics is 
observed, leading to an increase of the activation energy for the lithium diffusion by a factor 
of 1.7. Furthermore, it has turned out to be the immobilization of the oxygen particles 
that contributes most to this behavior. Analysis of suitable three-time-correlation functions 
has shown that forward-backward correlations of the lithium dynamics is much more pro- 
nounced. Thus, the minima of the potential landscape inside the fixed network structure 
act on the ions like an effective harmonic potential for significantly larger displacements of 
the ions than this would be the case inside a network carrying out its natural dynamics. 
Furthermore, by immobilizing the network species the lithium dynamics gets more hetero- 
geneous. Recently, also the effect of immobilization of a few ions on the the dynamics of the 
remaining ions has been studied jl^ . 

Another important ingredient for the ion dynamics through the glass matrix concerns the 
availability of a vacant site in the neighbourhood of an ion. This problem gains particular 
importance if the number density of ions is close the number density of sites. From molecular 
dynamics simulations by Lammert et ai, Habasaki et al. and Vogel it has turned out that 
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in the lithium sihcate system investigated in this work the number density of sites is only 



few percent larger than the one of ions 0, 14 1 . This result is consistent with the previous 
observation of ion channels Thus, for most ions there is no adjacent vacant site to 

which it might hop. These ions have to wait for an adjacent ion to perform a hop leaving 
a vacant side before they can perform a hop on their own. Thus, a macroscopic charge 
transport requires some cooperative hopping mechanism. Such correlated subsequent hops 
have been observed by Cormack et al. simulating a sodium silicate system . For lithium 
silicate correlated subsequent hops have been found by Habasaki et al., especially for ions 
with high mobility . They traced this behavior back to the prevention of a 

correlated backward hop due to the re-occupation of the initial site. 

The goal of this work is to identify ion hops and sample the relevant dynamic information. 
Our analysis goes beyond the observation of individual hopping scenarios. By calculating av- 
erages over many different hopping events (of the order 10^) we are able to make statistically 
relevant statements about the nature of the dynamics. 

We have organized this paper as follows. In sec. |n] we describe the technical aspects of 
the simulation and present the details of our analysis. The results are shown in sec. IIIII 
Finally, conclusions are drawn and discussed in sec. IIVI 

II. TECHNICAL ASPECTS 



A. Simulation 

The potential energy of the lithium silicate system is chosen to be the sum of a Bucking- 
ham and a Coulomb pair potential: 

?7.,(r,,) = -i-^ - ^ + A,, exp(-5,,r.,). (2.1) 

The ind.ce. , and , denote the =pec.e. litMun, oxygen, and .licon^espectively. The 



potential parameters are listed in previous wor 



on ab initio calculations by Habasaki et qL_ 
2.16e of the MOLDY software package 



c by 



22,122 



Banhatti and Heuer j2j|. They are based 



2J]. Applying a slightly modified version 



25|, we have generated molecular dynamics (MD) 
trajectories in the NVT ensemble, i.e. at constant volume. The length of the elementary 
time step has been chosen to be 2 fs and periodic boundary conditions have been used. We 
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have simulated a system of = 1152 particles at a density p = 2.34 gcm~^ taken from 
experimental room temperature data f26[]. If the temperature is not explicitly noted the 
presented data refers to T = 980 K. 

We have also modified the MOLDY software to compute only the dynamics of a chosen 
subensemble of particles and keep the positions of the other particles constant. Doing so, 
we are able to take a configuration of our system at an arbitrarily chosen time and continue 
the dynamics only for the ions so that the network is artificially immobilized. 



B. Definition of lithium hops 

Qualitatively, a hop is identified as a relatively large displacement during a relatively 
short time. Thus, we take a minimum displacement Smin occurring during a time interval 2r 
as the criterion for an ion hop. As we are interested in hops between next neighbor positions 
the value of Smin ought to be somewhat larger than the first minimum distance rmin ~ 1.5 A 
of the self part of the van Hove correlation function and somewhat less than the average 
hopping distance do ~ 2.6 A The average time an ion resides outside a hthium site has 
been found to be 500 fs Q]. The total time of the hop, also including the time to leave 
the initial site and to enter its new site is somewhat larger. Here we choose r = 1 ps. The 
actual choice is uncritical (see below). 

Let us assume to to be the time when a hopping ion is crossing the saddle point between 
two adjacent sites. To identify the position r^i of this saddle as well as the position rjnit of 
the initial lithium site and Tend of the lithium site after the hop, it has turned out to be very 
useful to average the position of the ion at times to, to — t, to + t, respectively, over a short 
time. Since a priori to is not known we define for these averages: 

"si 

rsi(t) = ^^^^^ r(t + «tsiice), (2.2) 



«=-nsl 

"init 



rinit(i) = 7^ — r(t - r + itsiice), (2.3) 

2ninit + 1 . ^ 

«=-"init 
-j^ "end 

rcnd(t) = r(t + r + itsiice). (2.4) 

2nend + 1 . ^ 



^ ^cnd 



tsiice = 20 fs is the time between subsequent system configurations, stored during the MD 
run. Using Usi = 5 and riinit = riend = 25 the position of the saddle point is averaged over 
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200 fs, and the positions of the hthium sites are averaged over 1 ps. 

How to find an appropriate choice for to? For this purpose we evaluate eqs. ()2.2|) to ()2.4|) 
for all available times t and for every lithium ion. For an ion hop one expects: 

|-'"end J^init I ^ "^miri' l^'^) 

For a single hopping event one typically finds a set of subsequent points in time t which do 
all fulfil the condition ()2.5|1 . We identify the exact time to of the hop as the time t when 
the hopping lithium ion is half-way between the initial site and the new site. Thus, the best 
value to representing this particular ion hop is chosen by the criterion: 

IkendW - rsi(t)| - |rsi(t) - ri„it(t)|| = min. (2.6) 

The average hopping distance (|rend — i"si|) is closest to do with the choice Smin = 2 A. 
Why is r = 1 ps a good choice? An appropriate criterion for the choice of r is that the 
number of identified lithium hops is a maximum. It turns out that the number of identified 
jumps is basically constant for 0.7 ps < r < 1.3 ps. Analyzing a simulation run of 20 ns 
length at T = 980 K with this method we find more than 130000 ionic hops. 



C. Coordinate Frame 

We wish to study the densities of the ions and the network species (see section Hi D|) in the 
temporal and spatial surrounding of each identified lithium hop. To average over all hops 
we need to define frames of reference that provide a certain sensitivity to the characteristic 
changes we are interested in. First, we define a time scale relative to to: 

ts = t-to. (2.7) 

As we are interested in dynamic events near the saddle point and the initial site we introduce 
difference vectors relative to the position of the saddle point and the initial site, respectively, 
i.e. 

rs = r - rsi (2.8) 

and 

ri = r - Finit. (2.9) 
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Obviously, one could also define a position vector fe relative to the position of the ion site 
after the hop, but due to time reversal symmetry, this would not lead to any new information. 

It seems to be most instructive to describe the position vectors rs and rj in polar coor- 
dinates: 



rs 



Os 
\<PsJ 



(2.10) 



Here, the angle between the hopping vector r^nd — ^mit and the particular position vector is 
denoted by 9s and 9i, respectively, (ps and (pi describe the angle in the plane perpendicular 
to the hopping vector. For reasons of isotropy no dependence on 0g or (pi is expected. The 
definitions are visualized in fig. ^ 



D. Definition of the investigated density functions 

Let {pi (rs, Os,ts)) denote the number density of particles at time ts belonging to species 
i inside a small volume element characterized by rs and ^s and averaged over all identified 
lithium hops. This quantity has been investigated for the hopping lithium ion (HLi) and 
the non-bridging oxygens (NBO). Actually, we also analyzed the bridging oxygens and the 
adjacent lithium ions, but these results will not be shown explicitly in this paper. By 
relating the values of (pi (rs, 6's, ts)) to the total number density pu of lithium or po of 
oxygen particles, respectively, one ends up with: 

PHLi(rs,6's,ts) = {pHu{rs,Os,ts)) /pu, 

(2.11) 

PNBo{rs,Os,ts) = {pNBo{rs,Os,ts)) /po ■ 



III. RESULTS 



A. Lithium dynamics 

The quantity PHLi('^S) ^s? ^s) is plotted in fig. |21for five different times —2 ps < ts ^ 0. By 
definition the contributing ions perform a hop from the left to the right. The time evolution 
is as follows. 1.) At ts = ^2 ps Puu shows a pronounced maximum at the initial site around 
rs ~ 1.2 A and 6s = it. From there an area of high density of hopping lithium spans the 
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saddle along the axis = and reaches into the new site up to rs ^ 1 A. This reflects 
the presence of lithium ions which are about to hop from the site at Tend to ri^^ and, of 
course, will jump back at ts = 0. Thus this intensity reflects back-and-forth dynamics. 2.) 
At ts = —500 fs the area where the density of hopping lithium is flnite cumulates around 
the initial site. The maximum of Phli has moved into the direction of the saddle and is now 
located at rs ~ 1 A and = ^r. 3.) At ts = —120 fs Phli shows a broad maximum along 
the direction of the hopping vector between the initial site and the saddle point with center 
at rs ~ 0.7 A and = ^r. The density is strongly concentrated around this maximum. 4.) 
At ts = ~40 fs the maximum concentration is close to rs ~ 0.4 A and ^s = ^- 5.) Finally, 
at ts = the density is radially distributed around the saddle point. It is maximal within 
a radius rs ~ 0.3 A and decays rapidly for larger values of rs. Due to the fact that the 
definition of rs\ involves an average over several time steps, finite values of rs are possible 
for to = 0. 

Of course, the cumulation of the density of hopping lithium while ts approaches zero 
is due to the fact, that by definition an ion is crossing the saddle close to ts = 0. The 
broadening at earlier times shows that different ions behave quite differently in the time 
interval just before the hop. 

Emphasizing the effect of highly occupied lithium sites, it has turned out that at least 
for a short time during the shown time interval —2 ps < ts < the site close to r^nd is 
occupied by another ion for about 83 % of all identified lithium hops. This does not only 
imply a high probability to find an ion at a site almost immediately before another ion 
will perform a hop to this site. It also means a high probability for a site to be occupied 
again by an ion just after another ion has hopped out of it. We found the time scale of this 
re-occupation to increase slightly with decreasing temperature, but even at the significantly 
lower temperature T = 640 K more then 65 % of the initial sites have been re-occupied by 
another ion at ts = 2 ps (not shown). 

B. Oxygen dynamics around the saddle 

Before we focus on the oxygen dynamics during an ion hop, we show the partial pair 
correlation functions 5'LiNBo('") for lithium and non-bridging oxygen (fig. 5'LiNBo('") is 
different from zero for distances r > 1.6 A. A pronounced maximum is found at r ^ 2.0 A 
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representing the first shell of non-bridging oxygen around lithium. The maxima for the 
second and higher order shells of non-bridging oxygen are less pronounced. 

We can now compare this distance with characteristic distances rs of Pnbo('^S) ^s? ^s)- 
The time evolution of Pnbo is shown in fig. 0] 1.) At ts = —2 ps the density of non- 
bridging oxygen indicates two minima around the centers of the two lithium sites. The 
minimum around the initial site seems to be a little more pronounced. This matches the 
slightly but clearly higher lithium density at this site. Pnbo shows a broad maximum around 
rg ^ 1.8 A and ~ i^jl. The shortest distance to the saddle where Pnbo reaches zero is 
at rs ~ 0.6 A. 2.) While time approaches zero the minima in non-bridging oxygen density 
get more and more pronounced. At the same time the density spike for ~ separating 
the two minima vanishes. The maximum of Pnbo broadens in 6's direction and narrows 
in rg direction. Its peak increases and moves to slightly higher values of rg. At tg = 
nearly everywhere inside an approximately spherical volume with radius rg ^ 1.5 A around 
the saddle the density of non-bridging oxygen has vanished. The maximum is shifted to 
rg ^ 1.9 A. A similar behavior is observed for the bridging oxygens. 

The time evolution of PnboI'^s, 6'g, tg) gives evidence that during an ion hop the oxygen 
particles around the saddle point significantly retreat. This mechanism can be compared 
to a sliding door which has to be opened to let the ion pass. To investigate this effect in 
more detail, we show another plot of Pnbo? in which we concentrate on a thin slice of 0.5 A 
thickness perpendicular to the hopping vector around rg = 0; see again fig. [H The results 
are shown in fig. We have also included the corresponding data for a system with an 
immobilized network at the same temperature T = 980 K. 

As the most pronounced effect oxygen density for rg < 1.5 A is dramatically reduced when 
tg approaches zero. In case of the immobilized network the density of NBOs is also negligible 
for rg < 1.5 A. Fig. El again shows the sliding door mechanism as already mentioned in the 
discussion of fig. 01 It seems to be related to a critical minimum distance between oxygen 
and the saddle at the time it is crossed by the hopping ion. This is rg ^ 1.5 A for non- 
bridging oxygen. Very close to this value the respective pair correlation functions decays to 
zero (see fig. EI)- The remaining small differences do not necessarily mean that during a hop 
the distance between lithium and non-bridging oxygen falls below 1.6 A, because the ion 
does not have to permeate the plane around the saddle strictly at rg = 0. 

To investigate the sliding door mechanism somewhat closer we analyze the distribution of 
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non-briding oxygen atoms at ts = for different subensembles. Subensemble Ai is defined as 
the set of oxygens with distances rs < 1.5 A from a saddle at time ts = —1 ps. Subensemble 
A2 is defined via the same criterion, but at time = —2 ps. The oxygen-distribution of 
subensemble A2 at time ts = ~1 PS is distributed around rg ^ 1.5 A and is therefore very 
different to subensemble Ai at that time (rs < 1.5 A); see fig. El This difference reflects 
the typical fluctuations of the oxygen atoms during 1 ps without an ion hop. Now the 
interesting question emerges whether the oxygen distribution at ts = ps is different for 
both subensembles. As also shown in fig. IHl both distributions are identical. Comparison 
with Pnbo('"s,^s) for all non-bridging oxygens at time ts = shows that this distribution 
is very similar to the overall distribution at ts = after appropriate scaling to the number 
of contributing particles. Thus on average the position of oxygens close to the saddle at 
time ts = 0, i.e. the time of the lithium hop, does not depend on its position in the past. 
Furthermore it turns out that the oxygen distribution Pnbo(^S) = 0) only weakly depends 
on temperature, as shown from the comparison of T = 980 K K with T = 640 K. The 
implications of these observations will be discussed below. 

C. Oxygen dynamics around the initial site 

In fig. Owe show the time evolution of the relative density Pnbo('^i,^s) of non-bridging 
oxygen inside a hemisphere around rinit covering angles tt /2 < 61 < n. At time ts = ~1 ps 
the first maximum for Pnbo(^I)^s) is found at rj ^ 2.1 A. The density decays to zero 
at ri ^ 1.4 A. When time tg approaches zero, i.e. the ion leaves the initial site, the 
corresponding peak gets broader and slightly shifts to larger distances. This trend stops 
around ts ~ 120 fs. For larger times the density tends to approach the initial density 
before the jump. The reason is that very likely the initial site will be populated by another 
ion. In summary, two major effects result after vacating an ion site. First, the network 
structures somewhat loosens as reflected by the broadening of the distribution and, second, 
the surrounding oxygen atoms slightly increase their distance from this site. 
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IV. CONCLUSIONS 



In this paper we have introduced a useful tool to identify the time and the position of 
ion hops. This enables us to study the temporal and spatial evolution of the relative density 
distribution of hopping ions and adjacent oxygen atoms. This leads to an instructive picture 
of the microscopic ion hopping dynamics. It turns out that there is a high probability for 
a site to be re-occupied by an ion almost immediately after it has been left. Even for the 
lowest investigated temperature T = 640 K this re-occupation takes place on a time scale 
of picoseconds. This behavior is similar to results obtained for a sodium silicate system by 
Cormack et al. [16]. It also matches the perception of Habasaki et al. Q,^, 19, 2ol |. 

We have shown that oxygen atoms near the saddle point of an ion hop retreat when the 
hop occurs. The oxygen motion reminds of a sliding door that opens to let the ion pass. 
Here a minimum distance to the saddle is kept by the oxygen which is very close to the 
distance where the corresponding pair correlation function with lithium vanishes. 

This cooperative dynamics between lithium and those oxygen atoms close to the saddle 
can be easily represented in an effective 2D-energy landscape of the type as sketched in 
fig- El Of course, in reality the energy landscape is 3 N- dimensional (N: particle number). 
The common reaction coordinate contains the lithium hop as well as the back-and-forth 
dynamics of the oxygen. Of course, for different saddles and different oxygen atoms the 
specific energy landscape may somewhat look different. At ts = one may define an 
effective oxygen potential Vi{Ro)- The index i counts the different effective potentials, 
experienced by the different oxygen atoms. Now two length scales can be determined, ^hom 
denotes the homogeneous length scale describing the range of possible fluctuations in the 
potential Vi{Ro). Close to the minima one expects in harmonic approximation ^hom oc 
T^/^. In contrast, ^het accounts for the heterogeneity when comparing hopping processes 
at different locations or different oxygen atoms. More specifically, it denotes the range of 
positions for which the different potentials Vi{Ro) display their respective minima. The two 
extreme limits ^het ^ ^hom and vice versa are sketched in fig. IHl To elucidate the ratio of 
both length scales in the present case we would like to repeat two important observations, 
presented above. First, for oxygen atoms close to the saddle the position at ts = does not 
depend on its position 1 ps before the hop as derived from the comparison of the behavior 
of two different subensembles. Second, the oxygen distribution at ts = only shows an 
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extremely weak temperature dependence. Both observations are compatible with the limit 
^het ^ ^hom- In this casc the position of any oxygen at ts = is very well characterized 
by its effective potential Vi{Ro) and does (in the extreme limit ^hom 0) not depend 
on its previous positions. Furthermore, since the effective energy landscape is to a large 
extent temperature independent one would not expect a major temperature dependence of 
the oxygen distribution, in agreement with observation. In any event, one may estimate 
^hom < 0.5 A, i.e. the width of PNBo('^s,is = 0). 

These results also reveal possible mechanisms for the dramatic slowing down of the ion 
dynamics upon immobilization of the network [3]. From the fig. EJone can see that many 
oxygen atoms have to move as much as 0.5 A to arrive at their optimum position at ts = 0. 
In case of an immobilized network this motion is not possible. This means that many 
oxygen atoms may be much further away from their optimum position during the transition 
as compared to ^hom- This automatically implies that the effective barrier for the combined 
lithium-oxygen transition will be much higher. This effect is sketched in fig. |H1 where the 
suppression of the oxygen dynamics will lead to a significantly higher barrier of the transition 
in probably most of the hopping events. Of course, if the oxygen distance to a saddle is 
less than 1.5 A the resulting saddle is so high that effectively no lithium hopping is possible 
along this path. One can estimate from the oxygen distribution in fig. Elat ts = —2 ps that 
this is the case for roughly 10% of all saddles. In any event, based on these considerations 
the dramatic increase of the activation energy upon immobilization can be understood. 

As a gedankenexperiment one may wonder for which (unphysical) oxygen mass one has 
the transition between the mobile and the immobile limit. Based on the above interpretation 
of the correlated hopping motion it is required that the oxygen can follow the lithium ion 
during the time of a hop, which takes of the order of 500 fs. Thus the velocity of the oxygen 
atom should be decorrelated on time scales faster than 500 fs. For this purpose we analyze 
the quantity w{t), which is the time derivative of the mean square displacement and twice 
the time integral of the velocity autocorrelation function 0,0- It is plotted in fig. Ell The 
short-time peak of w{t) is a measure for the typical velocity decorrelation time. It is of the 
order of 30 fs for oxygen and thus much faster than the typical jump times of the lithium 
ions. Thus one can estimate that by increasing the oxygen mass by (500/30)^ the slowing 
down of the lithium dynamics should start. 

Actually, a similar estimation has been performed by Sunyer et al. for MD simulations of 



11 



a sodium silicate system [5[ . They explicitly observed the slowing down in dependence of the 
(artificial) mass of the network atoms. To rationalize the slowing down they have compared 
the fourier transform of the velocity correlation function of the sodium atoms with that of 
the oxygen atoms. This quantity expresses the vibronic density of states. They found that 
an overlap is important for the sodium dynamics. As soon as the oxygen atoms become too 
slow and the overlap decreases the sodium diffusion constant starts to decrease, too. This 
happens if the oxygen mass is scaled by a factor of the order of 100-1000. This estimation 
is close to the present case. 

A further conclusion about the interrelation of lithium and network dynamics, as eluci- 
dated in this work, can be drawn from fig. ^| First, it is important to realize that around 
t ~ 1 ps there is a transition from localized lithium dynamics in a single site to hopping 
dynamics between adjacent sites. This has been explicitly shown in jSj. This time scale is 
also consistent with the present results because the hopping transition also roughly takes 
1 ps. One can directly observe that up to 1 ps the observable w{t) of the network atoms 
only weakly depends on temperature. This is in agreement with the expectation for the 
vibrational short-time limit { w{t) oc T^^^). Around 1 ps the temperature dependence dra- 
matically increases and is comparable with the temperature dependence of w{t) for the 
lithium ions. Since in this time regime the lithium ions start to hop this correlation may 
refiect the sliding door mechanism as well as the network relaxation at ion sites, discussed 
in this work. For these mechanism one indeed expects a strong relation between network 
and lithium dynamics. 

We just note in passing that the comparison of w{t) for different species also has impli- 
cations beyond the immediate topic of this work. In conductivity experiments one observes 
the sum of the transport of all charged constituents of the probe. Their contribution to cr{uj) 
is proportional to w{t ^ vr^/cj) P,!^]- Thus one can conclude that for t > 1 ps, i.e. in the 
hopping regime, the contribution is mainly due to lithium dynamics. 

In summary, we have elucidated in great detail the interrelation of lithium and network 
dynamics. It turns out that the reaction path for transition involves a significant oxygen 
dynamics. It can be characterized with a sliding door mechanism. There are strong indi- 
cations that this dynamics is crucial to reduce the barrier height as compared to the case 
where the oxygen dynamics is suppressed due to network immobilization. 
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Fig. U Definition of tlie coordinate frame used in this work. The thin shce around r^/ 
will be used further below. 

Fig. m Temporal and spatial evolution of the relative density Phli (^Sj ^s, ^s) of hopping 
lithium. 

Fig. 131 Pair correlation function gui'f') for lithium and non-bridging oxygen (z = NBO). 

Fig. m Temporal and spatial evolution evolution of the relative density Pnbo ('"Sj^Sj^s) 
of non-bridging oxygen adjacent to hopping lithium. 

Fig. m Temporal and spatial evolution of the relative density Pnbo (^S; ''^s) of non- 
bridging oxygen in a 0.5 A thin slice around the plane 6^ = n /2. 

Fig. IHl Relative density Pnbo(^S7^s) of non-bridging oxygen in a 0.5 A thin slice 
around the plane 6's = 7r/2 for different subensembles Ai and A2 compared 
to Pnbo (^Sj^s) for all non-bridging oxygens in this slice. Particles belonging to 
subensemble Ai are located at < 1.5 A at ts = ~1 PS- Particles belonging to 
subensemble A2 are located at rs < 1.5 A at ts = —2 ps. 

Fig. El Temporal and spatial evolution of the relative density PNBo(n5^s) of non- 
bridging oxygen in a hemisphere around Finit covering angles ti /2 < 9i< tx. 

Fig. El Schematic 2D energy landscape, reflecting the cooperative motion of the lithium 
and an oxygen atom. The thick line describes the reaction coordinate. The 
straight line for constant tq describes the transition path without oxygen dy- 
namics. In contrast, along the straight line for constant ru one can read off the 
effective potential Vi(i?o)for the oxygen atom at ts = 0. 

Fig. ini Sketch of the two extreme scenario of the effective potentials Vi{Ro), experienced 
by the oxygen atoms at time ts = 0. (a) Limit of strong curvature; (b) Limit of 
weak curvature. 

Fig. 1101 Time derivative w{t) of the mean square displacement for all particle species in 
the system at three different temperatures. 
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